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Abstract 

A rigorous demonstration that given appropriate data on our past light cone leads to the 
determination of the metric functions and all their time derivatives on our past light cone is 
presented, thus showing how to evolve the solution we obtain from data on the light cone off it in 
a well-defined and straightforward way. It also automatically gives a procedure for constructing 
the solution for all spherically symmetric, inhomogeneous cosmological Lemaitre-Tolman-Bondi 
models in observational coordinates as a Taylor series in time of however many terms we need. 
Our procedure takes into account the essential data giving the maximum of the observer area 
(angular-diameter) distance, and the redshift Zmax at which that occurs. This enables the 
determination of the vacuum-energy density which would otherwise remain undetermined. 

PACS numbers: 98.80.-k, 98.80.Es, 98.80.Jk, 95.36.+X 
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I. INTRODUCTION 



Since the pioneer work of Hubble in the twenties, it has been well known that the time 
span of cosmological observations is very small compared to the age of our Universe. Con- 
sequently, all we can aim to obtain on a cosmological time scale is data on one single light 
cone (our past light cone), which, for all practical purposes, can be considered the same as 
Bubble's. That means that we do not have the time variation of any cosmological data to 
start with. For this reason, observational cosmological modeling in which one wants to find 
the metric based only on observed quantities has been very much restricted to modeling our 
past light cone. That has been a handicap researchers in this field have been trying to over- 
come for some time - with very limited success - and a key issue for the future development 
of observational cosmology. 

Given this lack of knowledge of the time variations of cosmological data and assuming that 
the Universe is spherically symmetric around us, the only way of finding the time variation 
of a cosmological quantity is to derive a Taylor series. But that demands knowledge of the 
quantity itself and all its time derivatives on our past light cone, if we use observational 
coordinates. 

In this paper, we present a rigorous demonstration that, given appropriate data on our 
past light cone, we can find the metric functions and all their time derivatives on our 
past light cone, thus showing how to evolve the solution off it into the past or into the 
future. Besides being important in its own right, it also automatically gives us a procedure 
for constructing the solution as a Taylor series in time of however many terms we need. 
Moreover, our procedure fully takes into account an important piece of data, the maximum 
of the observer area distance, and the redshift at which it occurs. Without these extra 
observables, we do not have enough independent data to determine the model - in particular 
to determine the cosmological constant A. An intermediate and necessary step for doing 

n 

that is to show how an important result originally obtained by Bellaby [1] in the Lemaite- 
Tolman-Bondi (LTB) 3 + 1 formalism also follows from the observational cosmology field 
equations. 

In several papers Araiijo and Stoeger Araiijo et al. and {sl, Araiijo, Roveda and 



Stoeger j4j] demonstrated how to solve exactly the spherically symmetric (SS) Einstein 
field equations for dust in observational coordinates without assuming Friedmann-Lemaitre- 
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Rebert son- Walker (FLRW) and with cosmological data representing galaxy redshifts, ob- 
server area distances and galaxy number counts as functions of redshift. All these papers 
assumed that A = 0. In a recent paper, Araiijo et al. [5'] we demonstrated how this program 
may be carried out when A 7^ 0. These data are given, not on a space-like surface of constant 
time, but rather on our past light cone C~{pq)^ which is centered at our observational posi- 
tion pq "here and now" on our world line C. These results demonstrate how cosmologically 
relevant astronomical data can be used to determine the space-time structure of the universe 
- the cosmological model which best fits it. This has been the aim of a series of papers going 

n 

back to the Physics Reports paper by Ellis et al. [8]. The motivation and history of this " 
observational cosmology (OC) program" is summarized in Araiijo and Stoeger 2|]. 

In all these papers, we have shown that if the redshift, observer area-distance, and 
number-count data can be fit to FLRW functional forms (these are very special forms the 
data must take, if the universe is FLRW) then such data determines a bona fide FLRW 
universe - which is not a priori obvious. 

For that purpose, we have used an integration scheme that has been improved over the 
years, but still has had an element of uncertainty about it, in the sense that we did not have 
a precise mathematical way of determining the time evolution of one of the metric functions 
A{w,y) (see below) which is fundamental for the solution of the problem. That was due to 
the fact that apparently we have the gauge freedom of choosing A{w.,{]) in any way we like. 
However, our experience, acquired in treating all these cases, has shown that, unless A{w, 0) 
has a very specific functional dependency on w, the solution thus obtained would not satisfy 
the central conditions. Unless the central conditions are satisfied, we are not guaranteed 
that the null surface on which we assume we have the data is a past light cone of our world 
line. More specifically, we found that using the apparent gauge freedom of choosing A{w, 0) 
at will leads to a formal solution of the field equations which does not necessarily satisfy the 
necessary boundary conditions, namely, the central conditions. 

Thus, it has become clear that in fact we do not have a choice of choosing A{w, 0) since, in 
order to find a solution that also satisfies the central conditions we have, in a sense, to work 
backwards from the central conditions trying to figure out among all "possible" choices of 
A{w,Q) the one that leads to a solution satisfying the central conditions. We were successful 
in applying this procedure to various FLRW models but it was not satisfactory, because it 
is not an algorithmic procedure and depends on one's ability to spot the correct functional 
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dependency that leads to the desired solution. 

Our next step was to try to understand what causes this loss of gau^e freedom of choosing 
A{w, 0) and in Araiijo, Roveda and Stoeger and Araiijo et al. ^] we argued that the 
fulfillment of the following conditions: 

(1) A{wQ^y) is determined by the data and the central conditions; 

(2) The coordinate y is choosen to be a comoving radial coordinate; 

(3) The central conditions ([5]) 

removes the freedom of rescaling the time coordinate w and completely determines A{w,y). 
That is indeed a correct set of assumptions for proving that conjecture, but we were not able 
then to present a rigorous mathematical proof of the result. Here, as we shall demonstrate, 
we have been successful in proving that the more restricted set of assumptions, namely (1) 
and (3), leads to the desired result for all inhomogeneous spherically symmetric LTB models, 
thus filling this important gap in observational cosmology modeling. 

Recently, Hellaby and Alfadeel [gI, after rederiving the key OC equations in a somewhat 
different way and emphasizing the free functions which must be determined for their solution, 
developed an integration scheme for determining the time evolution of A{w^y) - as well as 
the other metric functions B{w,y) and C{w,y) - but only indirectly, by first determining 
the time function t{w,y) along all the matter world lines, and then finding A{w,y) from 
that A = tu, (see Eq. f lTSj) below). This is because they could not find an evolution equation 
for Aw, as they themselves emphasize. Here we demonstrate how to do that, which leads 
to a much more direct and streamlined integration scheme - without having to resort to 
integrating over the 3 + 1, {t,y) domain first. 

Furthermore, as mentioned earlier, we show how the data (especially the redshift and 
the observer-area - or angular-diameter - distance data) in conjunction with the central 
conditions eliminate the apparent gauge freedom and completely determine the solution 
both on our past light cone, and off it into the past. 

The primary aim of OC program is to strengthen the connections between astronomical 
observations and cosmological theory. We do this by allowing observational data to deter- 
mine the geometry of spacetime as much as possible, without relying on a priori assumptions 
more than is necessary or justified. Basically, we want to find out not only how far our ob- 
servable universe is from being isotropic and spatially homogeneous (that, is describable by 
an FLRW cosmological model) on various length scales, but also to give a dynamic account 
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of those deviations (Stoeger et al 




By using observational coordinates, we can thus formulate Einstein's equations in a way 
which reflects both the geodesic flow of the cosmological fluid and the null geometry of 
C (Po)) along which practically all of our information about the distant reaches of our 
universe comes to us - in photons. In this formulation the field equations split naturally 
into two sets, as can be easily seen: a set of equations which can be solved on C~{pq), that 
is on our past light cone, specified hy w = wq, where w is the observational time coordinate; 
and a second set which evolves these solutions off C~{po) to other light cones into the past 
or into the future. Solution to the first set is directly determined from the data, and those 
solutions constitute the "initial conditions" for the solution of the second set. 

In this paper, for completeness, we review some aspects of the problem of determining 
the solution of the exact spherically symmetric Einstein equations for dust in observational 
coordinates. We refer the reader to Ellis ei a/ j^, Kristian and Sachs [9i], Araiijo and Stoeger 
[2I and references therein for a complete account of the philosophy and the foundations of 
the OC approach leading to the integration of Einstein field equations in observational 
coordinates. 

In the next section we define observational coordinates, write the general spherically 
symmetric metric using them and present the very important central conditions for the 
metric variables. Section IIIII summarizes the basic observational parameters we shall be 
using and presents several key relationships among the metric variables. Section HVl presents 
the full set of field equations for the spherically symmetric case, with dust and with A 7^ 0. 
In Section |Vl we present a general integration scheme for all inhomogeneous spherically 
symmetric LTB models and in section |VT] we briefly discuss our conclusions. 



II. THE SPHERICALLY SYMMETRIC METRIC IN OBSERVATIONAL COOR- 
DINATES 

We are using observational coordinates (which were first suggested by Temple [l^). As 
described by Ellis el al [8] the observational coordinates = {w, y, 6, </)} are centered on 
the observer's world line C and defined in the following way: 

(i) w is constant on each past light cone along C, with u"'daW > along C, where is the 
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4- velocity of matter {u^-Ua = — 1). In other words, each w — constant specifies a past fight 
cone along C. Our past light cone is designated as w — wq. 



(ii) y is the null radial coordinate. It measures distance down the null geodesies - with affine 
parameter u - generating each past light cone centered on C. y = on C and dy/du > on 
each null cone - so that y increases as one moves down a past light cone away from C. 

(iii) 9 and are the latitude and longitude of observation, respectively - spherical coordi- 
nates based on a parallelly propagated orthonormal tetrad along C, and defined away from 
C by k'^daO = k'^dacj) = 0, where /c" is the past-directed wave vector of photons (/c"/Ca = 0). 

There are certain freedoms in the specification of these observational coordinates. In w there 
is the remaining freedom to specify w along our world line C. Once specified there it is fixed 
for all other world fines. There is considerable freedom in the choice of y - there are a large 
variety of possible choices for this coordinate - the affine parameter, z, the area distance 
C{w, y) itself. We normally choose y to be comoving with the fluid, that is u°'day = 0. Once 
we have made this choice, there is still a little bit of freedom left in y, which we shall use 
below. The remaining freedom in the 9 and 4> coordinates is a rigid rotation at one point 
on C. 

In observational coordinates the Spherically Symmetric metric takes the general form: 



where we assume that y is comoving with the fluid, so that the fluid 4- velocity is = A~^6!^. 

The remaining coordinate freedoms which preserves the observational form of the metric 
is a scaling of w and of y: 



The flrst, as we mentioned above, corresponds to a freedom to choose w as any time 
parameter we wish along C, along our world line at y — 0. This is usually effected by 
choosing A{w,0). As we have briefiy indicated already, this freedom disappears when we 
apply the central conditions (see below) . The second corresponds to the freedom to choose y 
as any null distance parameter on an initial light cone - typically our light cone at w — Wq. 



ds^ = -A{w, yfdw'^ + 2A{w, y)B{w, y)dwdy + C{w, yfdn'^, 



(1) 
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Then that choice is effectively dragged onto other hght cones by the ffuid ffow. y is comoving 
with the ffuid 4-velocity, as we have already indicated. We shall use this freedom to choose 
y by setting: 

MwQ^y) = B{wQ,y). (3) 

We should carefully note here that setting A{w, y) = B{w, y) off our past light cone w = wq 
is too restrictive. 

In general, these freedoms in w and y imply the metric scalings: 

A-.i=^A, 5^5 = ^5. (4) 
dw ay 

It is important to specify the central conditions for the metric variables A{w^y)^ B{w,y) 
and C{w,y) in Eq. ([T]) - that is, their proper behavior as they approach y = 0. These are: 

as y ^ : A{w, y) A{w, 0) 0, 
B{w,y)^B{w, 0)^0, 

C{w,y)^B{w,0)y = 0, (5) 
Cy{w,y) B{w,0). 

These important conditions insure that C, our world line, is regular - so that all functions 
on it our bounded, and that the spheres {w, y = constant) go smoothly to C as y ^ 0. They 
also insure that the null surfaces w = constant are past light cones of observers on C (See 



reference 



, especially section 3.2, p. 326, and Appendix A for details). 



III. THE BASIC OBSERVATIONAL QUANTITIES 

The basic observable quantities on C are the following: 

(i) Redshift. The redshift z at time Wq on C for a comoving source a null radial distance 
y down C~{po) is given by 

i+-4^i. (6) 

A{wo,y) 

This is just the observed redshift, which is directly determined by source spectra, once 
they are corrected for the Doppler shift due to local motions. It is consistent with and 
complements the first of central conditions in Eq. ([5]). 
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(ii) Observer Area Distance. The observer area distance, often written as tq, measured 
at time wo on C for a source at a null radial distance y is simply given by 



ro = C{wo,y), 



(7) 



provided the central condition ([5]), determining the relation between C{w,y) and B{w,y) 
for small values of y, holds. This quantity is also measurable as the luminosity distance di 
because of the reciprocity theorem of Etherington Uj (see also Ellis 121]), 



{l + zfC{wo,y). (8) 



(iii) The Maximum of Observer Area Distance. General 
maximum Cmax for a relatively small redshift Zmax (Hellaby 
and Araujo and Stoeger [3]). At Cmax, of course, we have 



y speaking, C{wo,y) reaches a 
1|; see also Ellis and Tivon [l^ 



dCjwo, z) ^ dC{wo,y) ^ ^ 
dz dy ' 

further conditioned by 

fCi^ < 0. (10) 
dz^ 

Furthermore, of course, as we shall review below, with the solution of the null Raychaudhuri 
equation (Eq. ( IMl) below), the data set will give us y = y{z), from which we shall be 
able to find ymax = Vmaxizmax)- These Cmax and Zmax data provide additional independent 
information about the cosmology. Without Cmax and Zmax we cannot constrain the value of 
A. 

(iv) Galaxy Number Counts. The number of galaxies counted by a central observer out 
to a null radial distance y is given by 

N{y)=A-K ij{wo,y)m B{wo,y)C{wo,y) dy, (11) 
Jo 

where fi is the mass-energy density and m is the average galaxy mass. Then the total energy 
density can be written as 

dz 1 

IJ,{wo,y) = mn{wo,y) = Mo{z) — — -, (12) 

dy B{wo,y) 
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where n{wo,y) is the number density of sources at {wo,y), and where 



^'"^7 In Tilt- 



Here dQ is the sohd angle over which sources are counted, and J is the completeness of the 
galaxy count, that is, the fraction of sources in the volume that are counted is J. The effects 
of dark matter in biasing the galactic distribution may be incorporated via m and/or J . In 
particular, strong biasing is needed if the number counts have a fractal behaviour on local 



scales (Humphreys et al 15(1 )• Iii order to effectively use number counts to constrain our 
cosmology, we shall also need an adequate model of galaxy evolution. We shall not discuss 
this important issue in this paper. But, fundamentally, it would give us an expression for 
m = m{z) in Eqs. (fT2l) and (fT3!) above. 

There are a number of other important quantities which we catalogue here for complete- 
ness and for later reference. 

First, there are the two fundamental four-vectors in the problem, the fluid four-velocity 
u"' and the null vector k"-, which points down the generators of past light cones. These are 
given in terms of the metric variables as 

^ A-^6\ , k"" = {AB)-^6%. (14) 

Then, the rate of expansion of the dust fluid is 3H = VaU'^, so that, from the metric (1) 
we have: 

where a "dot" indicates d/dw and a "prime" indicates d/dy, which will be used later. For 
the central observer H is precisely the Hubble expansion rate. In the homogeneous (FLRW) 
case, H is constant at each instant of time t. But in the general inhomogeneous case, H 
varies with radial distance from y = on t = Iq. From our central conditions above (3), we 
find that the central behavior of H is given by 

as.^O: i/(„.,)^^» = i/(„.0). (16) 

At any given instant w = Wq along y = 0, this expression is just the Hubble constant 
Hq = H(wo,0) = Aq^Bq^(B)o as measured by the central observer. In the above we have 
also written Aq = A{wo, 0) and Bq = B{wq, 0). 
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Finally, from the normalization condition for the fluid four- velocity, we can immediately 
see that it can be given (in covariant vector form) as the gradient of the proper time t along 
the matter world lines: Ua = —t,a- It is also given by ([T]) and (IT^ as 

Ua = 9ahU^ = -Aw^a + By^a- (17) 

Comparing these two forms implies 

dt = Adw - Bdy ^ A = ty, , B = -ty, (18) 

which shows that the surfaces of simultaneity for the observer are given in observational 
coordinates by Adw = Bdy. The integrability condition of Eq. (fT8l) is simply then 

A' + B = 0. (19) 

This turns out precisely to be the momentum conservation equation, which is a key 
equation in the system and essential to finding a solution. 

IV. THE SPHERICALLY SYMMETRIC FIELD EQUATIONS IN OBSERVA- 
TIONAL COORDINATES 

Using the fiuid-ray tetrad formulation of the Einstein's equations developed by Maartens 
{igI and Stoeger et al one obtains the Spherically Symmetric field equations in obser- 
vational coordinates with A 7^ (see Stoeger et al [2] for a detailed derivation). Besides the 
momentum conservation Eq. (fT9|) . they are as follows: 

A set of two very simple fluid-ray tetrad time-derivative equations: 



- = -3^f- + ^l (21) 



C\ 6 

where /im again is the relativistic mass-energy density of the dust, including dark matter, 
and 

is a quantity closely related to fimoiu) = f^m{wo,y) (see Eqj30l) below). 
Equations (!20l) and (|2T|) can be quickly integrated to give: 
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(22) 



u{w, y) = [..{y) + - j -_ = -__ + __ + -[-^), (23) 

where uJo{y) = uj{wo,y) and the last equahty in (l23l) follows from the definition of u given 
above. In deriving and solving these equations, and those below, we have used the typical 
A equation of state, pa = — /^a, where pa and fiA = are the pressure and the energy 
density due to the cosmological constant. Both ujq and /io are specified by data on our past 
light cone, as we shall show. fiA will eventually be determined from the measurement of 

Cfnax and ZjjiQ^x- 

The fluid-ray tetrad radial equations are: 



f = ?(;f + fj-i^V.; (24) 

(ujoiy) + if)C'{wo,y)^ = -i/i^o B{wo,y) C\wo,y) + (25) 
i = §^-L + ^)AB. (26) 



The remaining "independent" time-derivative equations given by the fluid-ray tetrad formu- 
lation are: 



^ ^^ + (^+^)a^; (27) 



C CA\2 
B B A , o 1 . o 

B-BA-^ - 2^" • (2^' 
From Eq. (!25|) we see that there is a naturally defined "potential" (see Stoeger et al [3]) 
depending only on the radial null coordinate y - since the left-hand-side depends only on y, 
the right-hand-side can only depend on y: 

C r" 

Fiy) + (29) 

Thus, from Eq. itself 

^o(y) = - 7^77^7^ V f fJ^moiy) B{wo,y) C^{wo,y) F{y) dy. (30) 
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Connected with this relationship is Eq. (l23l) . which we rewrite as 



CC A C'^ AB AB 

+ 



(31) 



CC 2BC^ 2C2 
where Cq = C{wo, y). 

We can now proceed to recover a simple but very important observational relationship 
which will enable us to determine /ia. We begin by differentiating Eq. fl29|) with respect to 
w. This gives 



C' = b[f-^-^^{AC-AC). 
Substituting for uj in Eq. fl26l) from Eq. fl23|) gives, using Eq. fl29l) 



(32) 



C' = 



2C2 ^ ACBC 2\BC) 2 



From equations fl5^ and (1551) we obtain 



(33) 



^2 

Eq. (!34l) can be rewritten as 



yl3 



(34) 



dw 



3 0. (35) 

Hence, the expression within the square brackets depends only on the radial null coordinate 
y and we define 



2Miy) 



^-CiF^-1)-^). 



From Eq. (136!) we find that 

C 
A 



± 



2M{y) , /iAC2 



C 



+ (F' - 1) + 



1/2 



Substitution of Eq. ([37D into Eq. ([29]) gives 



(36) 



(37) 



B A ^ 

Now, for y = ymax, C = 0. Therefore, 



1/2 



(38) 
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QMmax + f^ACl,, - SCmax = 0. (39) 

Eq. ( !39l) originally obtained by Hellaby ij] in the 3+1 framework has to be considered a 
fundamental relation in Observational Cosmology, since it enables, from C{wo, Zmax) and 
Zmax measurements, the determination of the unknown constant /iA (see below). 

Returning now to the main thread of our solution scheme, we substitute for C /B in Eq. 
(!23|l from Eq. (l29ll and use of equations (l30ll and (136!) gives 

M{y) = -L(y) + ^]ci{y) = I f\mM B{wo,y) C'{wo,y) F{y) dy. (40) 



From equations (ITTl) and (HOll we find that the mass parameter M{y) is related to number 
counts N{y) as follows: 

Miy) = ^ r mN\y)F{y)dy =^ T M\y)F{y)dy. (41) 
Svr JO Stt jo 

where, M{y) = mN{y) is the total mass summed over the whole sky by a central observer 
out to a null radial distance y. 



Stoeger et al 



and Maartens et al [18] have shown that equations (l29l) and ( 1311) can be 



transformed into equations for A and B, thus reducing the problem to determining C: 

C 

^ " [F2 - 1 + 2M/C + (/iA/3)C2]V2 (^2) 
^ " F ± [F2 - 1 + 2M/C + (/iA/3)C2] 1/2 • 



The LTB form of the exact solution (Lemaitre 19|, Tolman [20], Bondi 2l|; and cf. 
Humphreys ^2^ and references therein) is obtained by integration of (H2l) along the matter 
flow y =constant using flTSl) 

t - T{y) = I ^p2_i^ 2M/C + (/iA/3)C2]i/2 ' 
where T{y) is arbitrary, and we identify 

F^ = l-kf, k = 0,±l. (45) 

iere / = /(y) is a function commonly used in describing LTB models in the 3 + 1 coordinates 
23 1 . Hellaby and Alfadeel have made clear in their paper that they were unable to find 
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a time evolution equation for A{w,y) without resort to integrating over the 3 + 1, it,y) 
domain first. It is precisely Eq. (jH]) above that is used to initiate their indirect procedure 
to find A{w,y). In the next section we find a time evolution equation for {A{w,y) and a 
path towards full integration off our past light cone which avoids this detour. 



V. THE GENERAL SOLUTION - TIME EVOLUTION OFF OUR LIGHT CONE 

In this section we describe in detail the general integration procedure that is applicable to 
all inhomogeneous spherically symmetric universe models - that is the only constraint. We 
do not know whether the universe is homogeneous or not. But the data gives us redshifts z, 
observer area distances (angular-diameter distances) tq^z), "mass source densities" Mo{z), 
and the angular-distance maximum Cmaxiwojz) at Zmax- It is important to specify the 
latter, because, as we have already emphasized, without them, we do not have enough 
information to determine all the parameters of the space-time in the A 7^ case. For 
instance, although we can determine C{wo, z) with good precision (by obtaining luminosity 
distances and employing the reciprocity theorem, equation ([8])) out to relatively high 
redshifts, at present we do not yet have reliable data deep enough to determine Cmax and 
Zmax- But this has just recently become possible with precise space-telescope distance 
measurements for supernovae la. 



Mustapha, et al [2J] have shown that a A = LTB model can fit any reasonable set 
of redshifts, observer-area-distance and galaxy- number-counts (or equivalently our Mq{z)) 



data. However, as Hel 



also Krasihski, et al 



aby [1] first recognized, and as we have been insisting here (see 
25l |. and references therein), CmaxiwQ, z^ax) data place an additional 
constraint on solutions, which may require a nonzero A. This is very important to determine. 
Even allowing for large-scale inhomogeneities, is there definite evidence for nonzero vacuum 
energy, or some other form of dark energy? 

In pursuing the general integration with these data, we use the framework and the 
intermediate results we have presented in Section [IVl Obviously, one of the key steps 
we must take now is the determination of the "potential" F[y), given by Eq. (!29l) . This 
was done in a similar way for A = by Araiijo and Stoeger [2], as indicated above. This 
means we need to determine C'{wo,y) and C{wo,y), which we now write as C'q and Co, 



14 



respectively. We also need A{wo,y)- We remember, too, that at on w = wq we have chosen 
B{wo,y) = A{wo,y), which we have the freedom to do. 



Clearly, Cq can be determined from the tq^z) = C{wq,z) data, through fitting, along 
with the solution of the null Raychaudhuri Eq. to obtain z = z{y) (Stoeger et al 
A{wo,y), too, is obtained from redshift data along with this same z{y) result. We pause 
here to mention that this latter result together with Eq. ([6]) allows us to determine A[wo, y) 
and its value on C, A{wo, 0). The unknown constant A{wo, 0) in Eq. and in the solution 
of Eq. ( 12^ cancel when we equate the two results, allowing us obtain an explicit expression 
for A{wo,y). Setting y = 0, we then find A{wo,0) precisely. In a way Eq. provides us 
with a central condition ai w = Wq. 

Co is somewhat more difficult to determine. But the procedure is straight-forward. 

We determine Cq by solving Eq. ( l26l) for it on w = wq. Using equations ([3]) and (fT9l) . we 
can write this now as: 

Co{y) Ao{y)Co{y) 

But, from Eq. (123|) we can write u!o{y) in terms of Co{y), C'Q^y), and Co{y)- So Eq. (146|) 
becomes: 

This is a linear differential equation for Co(y), where from data we know everything on our 
past light cone, w = Wq, (once the null Raychaudhuri Eq. has been solved) except Co{y) 
itself and fi\, which is a constant that can be carried along and determined subsequently 
from C{wo, Zmax) and z^ax measurements (see below). Thus, we can easily solve Eq. fll7|) 
for Co{y), which will also depend on the unknown constant fi\. Its general solution is given 
by: 



where K is an integration constant to be determined by a boundary condition that, in our 
case, is a central condition for C{w, y). From Eq. ([5]) we can easily find that 
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asy^O: C{w,y) ^ B{w,0)y = 0. (49) 

Since from Eq. ([5]) Co(?/) — as y — we find that the constant K has to be set to zero 
if Eq. f H9|) is to be satisfied. This procedure enables us to determine F{y), that obviously 
also depends on /iA. Our next step is to insert this result along with N'{y) into equation 
( |4T1) to determine the mass function M{y). Next we evaluate the mass function M{y) at 
ymax and plug the result into Eq. (1391) which becomes an algebraic equation for With 
this determination of /ia, we know Co(?/) completely, and can now determine F{y) from Eq. 
f l29|) . From here on, we can now follow the solution oS w = Wq for all w as we demonstrate 
in detail below. 

We now turn to deriving an equation for A{wQ,y), which leads to finding subsequent 
equations for all higher derivatives of the metric variable on our past light cone. This leads 
to our key result. 

Equation (l28l) can be rewritten as 

Substitution of equation (HM into (15D1) yields 



dw\A 



d_ 

dw 



d_ 

dy 



(In A) 



d_ 

dy 



d_ 

dw 



(In A) 



d_(A^ 
dy\A^ 



2uj AB + AB 
2 



(51) 



where we used the fact that w and y are independent coordinates. 
On our past light cone {w = Wq) the above equation reads 



d 






dy 


_A{wo,y)_ 





2t^o(l/) 



A\wo,y). 



(52) 



Therefore, its general solution is 
Ao{y) = Aoiy) 



2^0 (y) + -^f^moiy) 



Ao'iy) dy + Ci 



(53) 



where we have written Aoi^y) and A^i^y) for A{wQ,y) and A{wo,y) respectively. Ci is an 
integration constant that one would expect, in principle, to be determined by a central 
condition for Ao{y), namely, Ao{0). (As mentioned above, we have already determined 
Aq{0), the first term of our Taylor series for A{w,y).) However, this procedure does not fix 
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Ci, since the central condition for A{w, y) only specifies that it must not go to zero as y goes 
to zero. Therefore, we have to investigate if use of the other central conditions provides a 
way of fixing Ci. 

It is important to note that it is this equation which enables us to stay within the OC 
formalism during the integration, providing a much more direct and streamlined alternative 
to the Hellaby and Alfadeel ^ scheme. 

We observe that ujQ{y) is fully determined at this stage since it can be obtained from the 
first equality in Eq. (HO!) as 



M{y) 



6 



(54) 



Hence, we have shown that the data on our past light cone determines Aoiji), except 
for the constant Ci which still remains undetermined at this stage. Since we know Ao(?/) 
from the data, equations (1271) and fl28|) evaluated on our past light cone become algebraic 
equations for Coiy) and Boiy)^ respectively 



" Co(y)Ao(y) 
<-^olZ/j - — -7-r\ — ^ 

A'o{y)My) 



My) + ^^mXy) 



Ao'{y)Co{y) 



My) 



Aoiy) 



Ao\y)Boiy) 



(55) 
(56) 



where in the later we have used Eq. f|T9l) . Note that Co(y) is completely determined at 
this stage except for its dependency on Ao{y) which carries its Ci dependency as we have 
explained above. However, Co{y) must satisfy the central condition for Co{w,y), that from 
Eq. © is 



as C{w,y)^B{w,0)y = 0, 



(57) 



which, applied on w = wq, fixes the constant Ci, and we are assured that Co{y) is completely 
determined by the data on our past light cone. 

Our next step is to differentiate Eq. ( l28l) with respect to w, that is. 



d^B Ab bA bA^ 



d_ 

dw 



2u; + -/im 



BA' 



(58) 



dw^ A A 

Substitution of Eq. (fT9l) and its second time derivative into the L.H.S. of Eq. (158|) gives 
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AA' 



A' 



BA BA^ 



A 

d_fA; 

dy\A^ 



A 



A^ 



d_ 

dw 



+ 



BA' 



BA BA^ 1 d 



1 

2' 



2uj + BA' 



(59) 



Evaluating Eq. fl59l) on our past light cone gives 



d_ 

dy 



My) 
Ao{y) 



i3,{y)A,{y) , Bo{y)Al{y) 



Al{y) 



+ 



1 [ d 

Ao{y)\dw 



My) 



2uj + 



BA' 



(60) 



1)0 



Therefore, its general solution is 



My) = -My)' 



y My)My) My)Al{y) 



Mij) 



My) ^9w 



2iu 



My) 

BA' 



dy + C2 



(61) 



where, using the same reasoning as above, C2 is an integration constant to be determined 
by the central condition for d'^C{wo,y), similar to Eq. (1571) . which is found as the solution 
of an algebraic equation once we differentiate Eq. (p71) with respect to w and substitute for 
the previously determined quantities. 

It is important to note that all quantities on the R.H.S. of the above equation are ob- 
tainable either directly from the data or from the algorithmic steps in the procedure we are 
describing here (Appendix). Therefore, we have shown that we can obtain A^ly) from the 
data. It is clear now that repetition of this procedure will give us all time derivatives of A, B 
and C on our past light cone, which means that A{w, y) B{w, y) and C{w, y) are completely 
determined by data on our past light cone, and calculable as Taylor series. 

It is clear from the above procedure that each step begins by finding the successive 
time derivatives of the metric function A{w,y) on our past light cone, d'!^A{wQ,y). For its 
complete determination one must apply the corresponding central condition, that is, one 
must specify 9^y4(t(7o, 0), which is done through other higher-derivative central conditions 
of the form given in Eq. (!57|) . Hence, it is shown that A{w, 0) is completely determined by 
the central conditions, thus proving the conjecture discussed in the introduction that the 
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fulfillment of assumptions (1) and (3) removes the freedom of rescaling the time coordinate 
w and completely determines A{w,y). 

VI. CONCLUSION 

We have summarized the essential details of previous work showing how to construct 
all spherically symmetric, inhomogeneous cosmological (LTB) models in observational 
coordinates from cosmological data on our past light cone, allowing for a nonzero cosmo- 
logical constant (vacuum energy). In doing so we provide a new rigorous demonstration 
of how such data fully determines the time evolution of all the metric components, 
and a Taylor series algorithm for determining those solutions. This enables us to 
move the solution we obtain from data on the light cone off it in a well-defined and 
straightforward way. It is essential for these to have data giving the maximum of the 
observer area (angular-diameter) distance, 6*0(^0, ^max), and the redshift Zmax at which 
that occurs. This enables the determination of the vacuum-energy density /^a, which 
would otherwise remain undetermined. That is, do we need A 7^ to adequately fit 
observational data, if we do not assume that the Universe is FLRW? Using this broader 
theoretical framework will enable us eventually to answer this important question, as well as 
to determine more securely how close or far the universe on large scales is from being FLRW. 

VII. APPENDIX 

The Algorithm: 

(i) Solve the null Raychaudhuri Eq. flM|) to obtain z = z{y). 

(ii) Determine Cg from the rQ{z) = C{wo, z) data, through fitting, along with the solution 
of the null Raychaudhuri equation. 

(iii) Determine Co by solving Eq. (|26ll for it on w = wq. Note that at this stage Co 
depends on /xa- 

(iv) Determine F{y) from Eq. (12^ which also depends on /ia. 

(v) Determine the mass function M[y) from equation (14T|) using F{y) along with N'[y) 
that is given from data. 
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(vi) Evaluate the mass function M{y) at Umax- 

(vii) Determine /^a from Eq. fl39|) . 

(viii) Knowing /ia, we determine both Co{y) and F{y) completely. 

(ix) Determine ujQ{y) from equation flS^ . From there on, we can follow the solution off 
w = Wq for all w. 

(x) Determine Ao^y) using the procedure described between equations (l50l) and (153|) 
applied to equation fl28l) . 

(xi) Determine both Co{y) and -Bo(l/) algebraically from equations fl27j) and (1281) . evalu- 
ated on our past light cone, respectively. 

(xii) Differentiate equations fl25]) (I27|) with respect to w. 

(xiii) Go back through steps (x) to (xii) to find at each run d^^A^wo, y) , d'^B{wQ, y) and 
dl^C{wo,y) respectively. 



[1] C.W. Hellaby, Mon. Not. R. Astron. Soc. 370, 239 (2006). 

[2] M.E. Araiijo and W.R. Stoeger, Phys. Rev. D 60, 104020 (1999). 

[3] M.E. Araujo, R.C. Arcuri, M.L. Bedran, L.R. de Freitas and W.R. Stoeger, Astrophys. J. 
549, 716 (2001). 

[4] M.E. Araujo, S.R.M.M. Roveda, and W.R. Stoeger, Astrophys. J. 560, 7 (2001). 
[5] M.E. Araujo, W.R. Stoeger, R.C. Arcuri and M.L. Bedran Phys. Rev. D 78, 063513 (2008). 
[6] C. Hellaby and A. H. A. Alfedeel, Phys. Rev. D 79, 043501 (2009). 
[7] W.R. Stoeger, G.F.R. Ellis, and S.D. Nel, Class. Quantum Grav. 9, 509 (1992). 
[8] G.F.R. Ellis, S.D. Nel, R. Maartens, W.R. Stoeger, and A.P. Whitman, Phys. Rep. 124, 315 
(1985). 

[9] J. Kristian and R.K. Sachs, Astrophys. J. 143, 379 (1966). 
[10] G. Temple, Proc. R. Soc. London A 168, 122 (1938). 
[11] I.M.H. Etherington, Philos. Mag. 15, 761 (1933). 

[12] G.F.R. Ellis, in General Relativity and Gravitation, edited by R.K. Sachs (Academic, New 

York,1971) p. 104. 
[13] G.F.R. Ellis and G. Tivon, Observatory 105, 189 (1985). 
[14] M.E. Araujo and W.R. Stoeger, Mon. Not. R. Astron. Soc. 394, 438 (2009). 



20 



[15] N. P. Humphreys, D. R. Matravers, and R. Maartens, Class. Quantum Grav. 15, 3041 (1998). 
[16] R. Maartens, Ph.D. thesis. University of Cape Town, 1980. 

[17] W. R. Stoeger, S.D. Nel, R. Maartens, and G.F.R. Ehis, Class. Quantum Grav. 9, 493 (1992). 
[18] R. Maartens, N. P. Humphreys, D. R. Matravers, and W. R. Stoeger, Class. Quantum Grav. 

13, 253 (1996); 13, 1689 (erratum) (1996). 
[19] G. Lemaitre, C. R. Acad. Sci. 196, 903 (1933); 196 1085 (1933). 
[20] R. C. Tolman, Proc. Nat. Acad. Sci. U.S.A. 20, 169 (1934). 
[21] H. Bondi, Mon. Not. R. Astron. Soc. 107, 410 (1947). 
[22] H. P. Humphreys, Ph.D. thesis. University of Portsmouth 1998. 
[23] W. B. Bonnor, Mon. Not. R. Astron. Soc. 167, 55 (1974). 

[24] N. Mustapha, C. Hehaby and G. F. R. Ehis, Mon. Not. Roy. Astron. Soc. 292, 817 (1997). 
[25] A. Krasihski, C. Hellaby, M-N. Celerier and K. Boleiko. [a?Xiv:0903.4070 vl [gr-qc] (2009). 



21 



